Deviation from quark number scaling of the anisotropy parameter V2 of pions, kaons, 
and protons in Au+Au collisions at ^/snn = 200 GeV 

A. Adare," S. Afanasiev,^^ C. Aidala,^^ N.N. Ajitanand,^^ Y. Akiba,50'5i H. Al-Bataineh ,^5 J. Alexander,^^ 
K. Aoki,32.50 Y. Aramaki,i° E.T. Atomssa,^^ R. Averbeck.^"^ T.C. Awes,-**^ B. Azmoun,^ V. Babintsev,22 
M. Bai,4 G. Baksay/8 L. Baksay/^ K.N. Barish,^^ B. Bassalleck,^* A.T. Basye,i S. Bathe,<^ V. Baublis,''^ 

C. Baumann,"" A. Bazilevsky,^ S. Belikov,^.^ R. Belmont,^! R. Bennett," A. Berdnikov,53 Y. Berdnikov,^^ 
A.A. Bickley," J.S. Bok,^^ K. Boyle," M.L. Brooks,^^^ H. Buesching,^ V. Buniazhnov,^^ G. Bunce,^'^^ S. Butsyk,^^ 

CM. Camacho,35 S. Campbell," C.-H. Chen," C.Y. Clii,!^ M. Cliiu,^ I.J. Choi,^^ R.K. Clioudliury,^ 
R Christiansen,^^ T. Chujo,60 P. Chung,56 0. Chvala,^ V. Cianciolo,'*^ Z. Citron," B.A. Cole,!^ M. Connors," 
P. Constantin,35 M. Csanad,!^ T. Csorgo,^'' T. Dahms," S. Dairaku,^^. so I. Danchev,^! K. Das,^^ A. Datta,^^ 
G. David,^ A. Denisov,^^ A. Deshpande,"' " E.J. Desniond,^ O. Dietzsch,^^ A. Dion," M. Donadelli,^"* 
O. Drapier,33 A. Drees," K.A. Drees,"* J.M. Durham," A. Durum,22 D. Dutta,^ S. Edwards,!^ Y.V. Efremenko,''^ 

F. Ehinghaus," T. Engelmore,*^ a. Enokizono,34 H. En'yo,^"'" S. Esumi,*^° B. Fadem,"** D.E. Fields,'^'' 
M. Finger,'^ M. Finger, Jr.,'^ F. Fleuret,^^ S.L. Fokin,^! Z. Fraenkel,'^^,^ j e. Frantz," A. Franz,^ A.D. Frawley,*^ 
K. Fujiwara,5° Y. Fukao,^^ T. Fusayasu,43 I. Garishvih,^^ A. Glenn," H. Gong," M. Gonin,33 Y. Goto,'5°'" 
R. Granier de Cassagnac, 33 N. Grau,i2 g y. Greene,^! M. Grosse Perdekamp,23." T. Gunji,i" H.-A. Gustafsson,37.0 
J.S. Haggerty,5 K.I. Hahn,i^ H. Hamagaki,!" J. Hamblen,^^ R. Han,^^ J. Hanks,*^ g p Hartoum,34 E. Haslum,37 
R. Hayano,i° X. He,20 M. Heffner,34 T.K. Hemmick," T. Hester,^ J.C. HiU,^^ M. Hohlmann,i« W. Holzmann,!^ 
K. Homma,2i B. Hong,30 T. Horaguchi,^! D. Hornback,^^ S. Huang,^! T. Ichihara,5°'" R. Ichimiya,^" J. Ide,"** 
Y. Ikeda,^" K. Imai,32.50 M. Inaba,^^ D. Isenhower,* M. Ishihara,50 T. Isobe,!"'^" M. Issah,^! A. Isupov,^^ 

D. Ivanischev,49 B.V. Jacak,"'[l| J. Jia,5^56 j^^i2 Johnson,^ K.S. Joo,^^ Jouan,47 D.S. Jumper,^ 
F. Kajihara,-^° S. Kametani,^^ N. Kamihara,^* J. Kamin,^'' J.H. Kang,^^ J. Kapustinsky,35 K. Karatsu,^^' 

D. Kawan,39." M. Kawashima,5°'^2 A.V. Kazantsev,3i T. Kempel,^^ A. Khanzadeev,'^^ K.M. Kijima,^! B.I. Kim,30 
D.H. Kim,'*2 D.J. Kim,28 E. Kim,^^ E.J. Kim,« S.H. Kim,^^ Y.J. Kim,23 E. Kinney," K. Kiriluk," A. Kiss,!^ 
E. Kistenev,^ L. Kochenda,^^ B. Komkov,''^ M. Konno,^" J. Koster,23 D. Kotchetkov,'''* A. Kozlov,^^ Krai," 
A. Kravitz,i2 G.J. Kunde,^^ K. Kurita,^*^' Kurosawa,5° Y. Kwon,^^ G.S. Kyle,^^ R. Lacey,^^ y.S. Lai,!^ 

r n T Qi^io 26 A T oKorloir 26 n A/T T oo 35 T T oo 17 XT T oo 55 T<- R T „ 30 T<- c T oo 30 iv/r T T oi + ^V, 35 i\/r A T T oHo 54 



A. Kravitz,i2 G.J. Kunde,^^ K. Kurita,^*^' Kurosawa,5° Y. Kwon,^^ G.S. Kyle,^^ R. Lacey,^^ y.S. Lai,!^ 
J.G. Lajoie,26 A. Lebedev,^^ D.M. Lee,^^ J. Lee," K. Lee,^^ K.B. Lee,^^ K.S. Lee,^^ M.J. Leitch,35 M.A.L. Leite,' 

E. Leitner,<^i B. Lenzi,^^ X. Li,^ P. Liebing,^! L.A. Linden Levy," T. Liska," A. Litvinenko,^'^ H. Liu,35.45 
M.X. Liu,3^ B. Love,6i R. Luechtenborg,'*^ D. Lynch,^ C.F. Maguire,^i Y.I. Makdisi,'* A. Malakhov,^^ M.D. Malik, 
V.I. Manko,3i E. Mannel,!^ Y. Mao,'*8.50 H. Masui,*^° F. Matathias,!^ M. McCumber," P.L. McGaughey,^^ 
N. Means," B. Meredith,23 Y. Miake,^^ A.C. Mignerey,38 p. Mikes,^'^^ K. Miki,5°'S0 Milov,^ M. Mishra,^ 
J.T. Mitchen,^ A.K. Mohanty,3 Y. Morino,*" A. Morreale,^ D.P. Morrison,^ T.V. Moukhanova,^^ J. Murata,^^'^^ 

S. Nagamiya,29 J.L. Nagle," M. Naglis,^^ j^j^ Nagy,!^ I. Nakagawa,^"'" Y. Nakamiya,^! T. Nakamura,^!^ 29 
K. Nakano,^°'5^ J. Newby,^'' M. Nguyen," R. Nouicer,^ A.S. Nyanin,^! E. O'Brien,^ S.X. Oda,i" C.A. Ogilvie,^^ 
M. Oka,s° K. Okada," Y. Onuki,^" A. Oskarsson,37 M. Ouchida,2i' ^o k. Ozawa,i" R. Pak,^ V. Pantuev,24." 
V. Papavassiliou,^^ I.H. Park," J. Park,^^ S.K. Park,30 W.J. Park,3" S.F. Pate,*^ H. Pei,26 J.-C. Peng,23 
H. Pereira,*^ V. Peresedov,2'' D.Yu. Peressounko,^^ C. Pinkenburg,^ R.P. Pisani,^ M. Proissl,^'' M.L. Purschke,^ 

A.K. Purwar,35 H. Qu,2" J. Rak,28 A. Rakotozafindrabe,33 I. Ravinovich,63 ^ p Read,46.58 K. Reygers,^^ 
V. Riabov,49 Y. Riabov,'^^ E. Richardson,38 D. Roach,^! G. Roche,^^ S.D. Rolnick,^ M. Rosati,26 C.A. Rosen," 
S.S.E. Rosendahl,37 p. Rosnet,^^ P. Rukoyatkin,27 p. Ruzicka,25 B. Sahlmueller,^" N. Saito,29 T. Sakaguchi,^ 
K. Sakashita,5"'59 y. Samsonov,^^ S. Sano,i°'62 T. Sato,^^ S. Sawada,29 K. Sedgwick,^ J. Seele," R. Seidl,23 

A.Yu. Semenov,26 R. Seto,^ D. Sharma,^^ shein,22 T.-A. Shibata,^". 59 K. Shigaki,2i M. Shimomura,60 
K. Shoji,32.50 p. Shukla,3 A. Sickles,^ C.L. Silva,^^ D. Silvermyr,^^ q Silvestre,!^ K.S. Sim,^^ B.K. Singh,2 
CP. Singh,2 V. Singh,2 M. Slunecka,^ R.A. Soltz,34 ^ Sondheim,35 S.P. Sorensen,^^ I.V. Sourikova,^ 
N.A. Sparks,! P.W. Stankus,'*^ E. Stenlund,^^ S.P. StoU,^ T. Sugitate,2i A. Sukhanov,^ J. Sziklai,^^ 
E.M. Takagui,54 A. Taketani,^^^ " ^ Tanabe,^" Y. Tanaka,43 Tanida,32.50.5i M.J. Tannenbaum,^ S. Tarafdar," 
A. Taranenko,56 p. Tarjan,!^ H. Themann," T.L. Thomas,*^ M. Togawa,32.50 A. Toia," L. Tomasek,25 
H. Torii,2i R.S. Towell,! I. Tserruya,63 Y. Tsuchimoto,2i C. Vale,5>26 h. Valle,^! H.W. van Hecke,^^ 
E. Vazquez-Zambrano,*2 a. Veicht,23 J. Vclkovska,''* R. Vertesi,*5,64 j^^ Vinogradov,3i M. Virius," 
V. Vrba,25 E. Vznuzdaev,*^ X.R. Wang,^^ D. Watanabe,2i K. Watanabe,6° Y. Watanabe,^^^ " F. Wei,^^ 
R. Wei,56 J. Wessels,''" S.N. White,^ D. Winter, ^2 J.p. Wood,* C.L. Woody,^ R.M. Wright, ^ M. Wysocki," 
W. Xie," Y.L. Yamaguchi,i° K. Yamaura,^^ R. Yang,23 A. Yanovich,22 J. Ying,20 S. Yokkaichi,^"'" Z. You,''^ 



2 



G.R. Young,46 J Younus,^^ I.E. Yushnianov,^! W.A. Zajc,!^ C. Zhang,^^ g zhou,^ and L. Zolin^^ 

(PHENIX Collaboration) 

'Abilene Christian University, Abilene, Texas 79699, USA 
^Department of Physics, Banaras Hindu University, Varanasi 221005, India 
^Bhahha Atomic Research Centre, Bombay 400 085, India 
Collider- Accelerator Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
^Physics Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
^University of California - Riverside, Riverside, California 92521, USA 
'''Charles University, Ovocny trh 5, Praha 1, 116 36, Prague, Czech Republic 
^ Chonbuk National University, Jeonju, 561-756, Korea 
''Science and Technology on Nuclear Data Laboratory, China Institute of Atomic Energy, Beijing 1024-13, P. R. China 
Center for Nuclear Study, Craduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

" University of Colorado, Boulder, Colorado 80309, USA 
'^Columbia University, New York, New York 10027 and Nevis Laboratories, Irvington, New York 10533, USA 
'"'Czech Technical University, Zikova 4, 166 36 Prague 6, Czech Republic 
'^Dapnia, CEA Saclay, F-91191, Cif-sur-Yvette, France 
'^Debrecen University, H-4OIO Debrecen, Egyetem ter 1, Hungary 
'^ELTE, Eotvos Lordnd University, H - 1117 Budapest, Pdzmdny P. s. 1/A, Hungary 
'^Ewha Womans University, Seoul 120-750, Korea 
'^Florida Institute of Technology, Melbourne, Florida 32901, USA 
'^Florida State University, Tallahassee, Florida 32306, USA 
^"Ceorgia State University, Atlanta, Georgia 30303, USA 
^'Hiroshima University, Kagamiyarna, Higashi- Hiroshima 739-8526, Japan 
^^IHEP Protvino, State Research Center of Russian Federation, Institute for High Energy Physics, Protvino, 142281, Russia 

University of Illinois at Urbana- Champaign, Urbana, Illinois 61801, USA 
^'^ Institute for Nuclear Research of the Russian Academy of Sciences, prospekt 60-letiya Oktyabrya 7a, Moscow 117312, Russia 
Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague 8, Czech Republic 

^^lowa State University, Ames, Iowa 50011, USA 
Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 
Helsinki Institute of Physics and University of Jyvaskyld, P.O.Box 35, FI-4OOI4 Jyvaskyld, Finland 
^^KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 
^"Korea University, Seoul, 136-701, Korea 
^'Russian Research Center "Kurchatov Institute", Moscow, 123098 Russia 
''^Kyoto University, Kyoto 606-8502, Japan 
Laboratoire Leprince-Ringuet, Ecole Polytechnigue, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
^'^ Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
^^Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
^^LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 
^''Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
University of Maryland, College Park, Maryland 20742, USA 
Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003-9337, USA 
'^"institut fur Kemphysik, University of Muenster, D-48149 Muenster, Germany 
'^'Muhlenberg College, Allentown, Pennsylvania 18104-5586, USA 
^^Myongji University, Yongin, Kyonggido 449-728, Korea 
Nagasaki Institute of Applied Science, Nagasaki-shi, Nagasaki 851-0193, Japan 
University of New Mexico, Albuquerque, New Mexico 87131, USA 
'^'"'New Mexico State University, Las Cruces, New Mexico 88003, USA 
^^Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
^''IPN-Orsay, Universite Paris Sud, GNRS-IN2P3, BPl, F-91406, Orsay, France 
Peking University, Beijing 100871, P. R. China 
■^■'PNPI, Petersburg Nuclear Physics Institute, Gatchina, Leningrad region, 188300, Russia 
''"RIKEN Nishina Center for Accelerator- Based Science, Wako, Saitama 351-0198, Japan 
^'RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
Physics Department, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima, Tokyo 171-8501, Japan 
Saint Petersburg State Polytechnic University, St. Petersburg, 195251 Russia 
^'^ Universidade de Sao Paulo, Instituto de Fisica, Caixa Postal 66318, Sao Paulo CEP05315-970, Brazil 

''Seoul National University, Seoul, Korea 
Chemistry Department, Stony Brook University, SUNY, Stony Brook, New York 11794-3400, USA 
^''Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, New York 11794-3400, USA 

^'^ University of Tennessee, Knoxville, Tennessee 37996, USA 
Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan 
^''institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
^' Vanderbilt University, Nashville, Tennessee 37235, USA 



3 



Waseda University, Advanced Research Institute for Science and 
Engineering, 17 Kikui-cho, Shinjuku-ku, Tokyo 162-0044, Japan 
^'^Weizmann Institute, Rehovot 76100, Israel 
Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Hungarian 
Academy of Sciences (Wigner RCP, RMKI) 11-1525 Budapest 114, POBox 49, Budapest, Hungary 
^^Yonsei University, IPAP, Seoul 120-749, Korea 
(Dated: February 26, 2013) 

Measurements of the anisotropy parameter V2 of identified hadrons (pions, kaons, and protons) as 
a function of centrality, transverse momentum pr, and transverse kinetic energy KEt at midrapidity 
{\r]\ < 0.35) in Au+Au collisions at -^/sjviv = 200 GeV are presented. Pions and protons are identified 
up to pt = 6 GeV/c, and kaons up to pT = 4 GeV/c, by combining information from time-of- flight 
and aerogel Cerenkov detectors in the PHENIX Experiment. The scaling of «2 with the number 
of valence quarks (n^) has been studied in different centrality bins as a function of transverse 
momentum and transverse kinetic energy. A deviation from previously observed quark-number 
scaUng is observed at large values of KET/wg in noncentral Au+Au collisions (20-60%), but this 
scaling remains valid in central collisions (0-10%). 

PACS numbers: 25.75.Dw, 25.75.Ld 



I. INTRODUCTION 

Measurements of the anisotropy parameter V2 (the sec- 
ond coefficient in the Fourier expansion of the hadron 
yields with respect to the reaction plane) have played 
a pivotal role in the discovery of the strongly coupled 
quark-gluon plasma (sQGP) at RHIC [H-Q. At low pr 
(< 2 GeV/c) the agreement between ideal hydrodynam- 
ics calculations and the data have led to the conclusion 
that a near-perfect fluid is created in heavy-ion collisions 
at RHIC Q • Recent theoretical efforts aiming to quan- 
tify the ratio of the shear viscosity to the entropy density 
77/s (see for example reviews in [zl-Q) have confirmed 
that in the sQGP fluid this ratio is close to a conjec- 
tured quantum limit [1Q|. The high pr (> 6 GeV/c) az- 
imuthal anisotropics [ll'-H^ have been attributed to the 
path-length dependence of energy loss in the medium and 
are used to constrain the theoretical descriptions of jet 
energy loss [H, [H. At intermediate pr (2-6 GeV/c), 
which is the focus of this paper, the identified hadron 
anisotropics have shown strong evidence for quark-like 
degrees of freedom and significant collectivity at the par- 
ton level. This is supported by the observation of scaling 
with the number of valence quarks in the hadron (n^ 
scaling) 

The scaling with number of valence quarks (n^) was 
seen as a confirmation of quark recombination as a novel 
particle-production mechanism that competes with frag- 
mentation in the intermediate-pT range. Recombination 
models [2^-|26j were developed to account for the unusu- 
ally large baryon-to-meson ratios (relative to V+P col- 
lisions) and nuclear- modification factors [itI . Wh . [28j , as 
well as the large elliptic flow at intermediate pt, with 
pronounced differences between baryons and mesons. 



'Deceased 

t PHENIX Spokesperson: | j acak® skipper . physics . sunysb. edu| 



In the models, the Uq scaling, which is manifested as 
yhadron^p^-^ « nqV2{pT / nq) , IS an approximate scaling 
that comes from the addition of the valence-quark mo- 
menta at hadronization, with the assumption that the 
collective flow develops at the partonic level. 

There are several theoretical considerations that sug- 
gest that the Uq scaling should be violated in certain 
conditions. For example, the inclusion of higher Fock 
states describing the contribution of sea quarks and glu- 
ons have been shown to affect the Uq scaling [29j. Simi- 
larly, models that consider recombination between "ther- 
mal" partons (soft partons thermalized in the medium) 
and "shower" partons (partons fragmented from jets) 
predict centrality-dependent deviations from riq scaling 
that are particle-species dependent [SOj. Understanding 
the limits of the recombination domain is important in 
relation to viscous hydrodynamics and the extraction of 
the shear viscosity over entropy density (77/s) from the 
data [sll - fsH , as well as for developing a unified approach 
in describing jet energy loss and high px V2 |34j3a |. 
Searches for deviations from Uq scaling are also impor- 
tant for the low-energy scan program at RHIC as they 
have been considered as a signature of the transition be- 
tween sQGP formation and a hadronic system. Recent 
considerations of baryon transport may complicate this 
picture [stJ, which further reinforces the need for a de- 
tailed understanding of this scaling at ^/snn = 200 GeV. 

The Uq scaling has been tested in certain centralities 
and Pt regions with identified particles [T6l - [2]| . However, 
the precision of experimental data on identified hadron 
V2 is in many cases limited in statistics and px reach, 
especially for baryon measurements at KEt /nq > 1 GeV 
(KEt = V Pt ~^ niQ — mo) where the riq scaling may start 
to break. Therefore, the detailed pt limits and centrality 
dependence of the Uq scaling have not been tested. 

This paper reports on high-statistics measurements of 
the second order Fourier coefficient V2 for identified pi- 
ons (7r+ + TT^), kaons {K^ + K^), and protons {p A-p), 
which extend to relatively high pT (up to 6 GeV/c for 
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pions and protons and 4 GeV/c for kaons). The data for 
different centrality events (0-10%, 10-20%, 20-40%, 40- 
60%, and combinations tfiereof) are analyzed separately 
and the riq scaling is examined as a function of centrality. 
Comparisons with published measurements of Kg and A 
from STAR collaboration [2l| are shown in the central- 
ities 0-10% and 10-40%. The experimental details are 
presented in Section UTJ the analysis methods are in Sec- 
tion |llTl the results and discussion are in Section HVl and 
Section fVl summarizes our findings. 

II. EXPERIMENTAL SETUP 

The PHENIX experiment is designed for the study of 
nuclear matter in extreme conditions through a variety of 
experimental observables. It comprises a tracking system 
optimized for the high-multiplicity environment of ultra- 
relativistic heavy ion collisions, a set of particle identifi- 
cation (PID) detectors, and a set of detectors aimed at 
determining the global properties of the collisions. 

Figure [1] shows a schematic diagram of the PHENIX 
detector. The upper part is a beam-axis view of the 
two central spectrometer arms (West and East), covering 
the pseudorapidity region of jryj < 0.35. Below that is a 
side view showing the two forward-rapidity muon arms 
(South and North) and the global detectors. A detailed 
description of the complete set of detectors can be found 
elsewhere [38{. 

The physics analysis presented here employed the 
tracking system [drift chamber (DC) and three layers of 
multiwire proportional chambers (PCI, PC2, and PC3)], 
the West arm time-of-flight detector (TOFw), the aero- 
gel Cerenkov counter (ACC) , the ring imaging Cerenkov 
counter (RICH), the beam-beam counters (BBC), the 
reaction-plane detector (RxNP), and the muon piston 
calorimeter (MPC). Below, we give a brief description 
of each of these detector sub-systems and their role in 
the present analysis. 

A. Global Detectors 

The BBC are located at ±144 cm from the nominal 
interaction point along the beam line in the pseudora- 
pidity region 3.0 < |ry| < 3.9. Each BBC comprises 64 
Cerenkov telescopes, arranged radially around the beam 
line. The BBCs provide the main collision trigger for 
the experiment and are used in the determination of the 
collision vertex position along the beam axis (z-vertex) 
and the centrality of the collisions. They also provide the 
start time for the time-of-flight measurement with timing 
resolution of ctbbc = 37 ps. 

The RxNP ^ was installed in PHENIX before the 
2007 data-taking period. It is located at ±38 cm from 
the nominal interaction point and has full azimuthal cov- 
erage. Each RxNP comprises two rings of plastic scintil- 
lator paddles, with each paddle subtending A0 = 7r/6. 




West Beam View East 




FIG. 1: (color online) The PHENIX detector configuration 
for RHIC 2007 data-taking period. The beam-beam counters 
(BBC) are labeled as BB. 

The inner and outer segments cover the pseudorapidity 
ranges 1.0 < \r]\ < 1.5 and 1.5 < \r]\ < 2.8, respectively. 
The RxNP is the main detector used in the event-plane 
determination for this analysis. The event-plane resolu- 
tion (Res(\E')) [131, which is used as a correction to the 
V2 measurement, is defined as 

Res(*) = (cos(2(«' - ^>Rp)) . (1) 

Here the bracket () indicates an average over all events, 
RP is true reaction plane (which is defined by the direc- 
tions of interaction parameter and beam), and is the 
event plane (which is measured by the detector event by 
event). A larger value of Res(^) corresponds to a better 
measurement of the event-plane. In a given event, the 
event-plane resolution depends on the charged- particle 
multiplicity and the size of the azimuthal anisotropy sig- 
nal; thus the resolution is centrality dependent. A reso- 
lution of up to 73% is achieved for midcentral events. 

The MPC are electromagnetic calorimeters situated at 
±223 cm from the nominal interaction point inside the 
cylindrical openings at the front of the muon magnet pis- 
tons [m and have 2it azimuthal acceptance. The pseu- 
dorapidity coverage is about 3.0 < rj < 3.8 for the north 
side and -3.7 < < -3.1 for the south side. The MPCs are 
comprised of 220 modules in the north piston hole and 
192 in the south with PbW04 crystals and Avalanche 
Photodiode readouts, and can detect both charged and 
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neutral particles. In this analysis, the MPCs were used 
for event-plane determination. Although the event-plane 
resolution (up to 50% in midcentral collisions) is lower 
than that achieved with the RxNP, the MPCs provide 
an important systematic check on the RxNP measure- 
ment due to their larger pseudorapidity separation from 
the central spectrometer and therefore smaller nonflow 
effects on the V2 measurement. 



B. Tracking and Particle Identification Detectors 

The charged particle momentum is reconstructed in 
the tracking system comprised of the DC located outside 
of an axially-symmetric magnetic field at a radial dis- 
tance between 2.0 m and 2.4 m followed by the PCI with 
pixel-pad readout. The pattern recognition in the DC is 
based on a combinatorial Hough transform in the track 
bend plane. A track model based on a field-integral look- 
up table determines the charged particle momentum, the 
path length to the TOFw and a projection of the track 
to the outer detectors. The momentum resolution in this 
data set was estimated to be Sp/p 1.3% ® 1.2% x p 
(GeV/c), where the first term represents multiple scat- 
tering up to the DC and the second term is due to the DC 
spatial resolution. The momentum resolution is worse 
than that of previous data set is due to the weaker mag- 
net configuration. 

The tracks are matched to hits registered in the sec- 
ond and third layers of the pad chambers, PC2 and PCS, 
which are located at radial distances of 4.19 m and 4.98 m 
from the interaction point. Thus, the contribution of 
tracks originating from decays and 7-conversions is re- 
duced. 

To improve the track purity further, we employ the 
RICH, which is a threshold gas Cerenkov detector located 
in the radial region 2.5 m < r < 4.1 m. The Cerenkov 
radiator gas (CO2) at atmospheric pressure has an index 
of refraction n — 1.000410 {jth — 35), which corresponds 
to a momentum threshold of 20 MeV/c for an electron 
and 4.65 GeV/c for a pion. The RICH provides a veto 
for the electrons and positrons, which are predominantly 
pairs resulting from 7-conversions and Dalitz decays. 

The primary PID used in this analysis is the TOFw, 
which is located at a radial distance of 4.81 m from the 
interaction point and covers the pseudorapidity range 
|7;| < 0.35 and Sc/} = 22° in azimuth. The TOFw was 
built using inultigap resistive plate chamber technology 
(MRPC) '4^ and installed in PHENIX before the 2007 
data taking period. The MRPCs have six gas gaps 
formed by layered glass plates with thickness of 550 /im, 
separated by 230 /im-thick monofilament fishing line. 
The MRPCs are positioned in a gas volume and oper- 
ated with a gas mixture of 95% R134a and 5% isobutane 
(C4H10), and bias voltage of 14 kV. The TOFw system 
is composed of 128 MRPC modules each of which has 
four signal strips of size 37 x 2.8 cm^ and separation of 



for hit positioning along the direction of the strip to be 
determined using the timing difference between the sig- 
nals with resolution of the order 1 cm. The other two hit 
coordinates are determined using the global position of 
the strips within PHENIX. The average time measured 
on both sides of the strips provides the stop time for the 
time-of-fiight measurement. The timing resolution of the 
BBC-TOFw system was determined by selecting charged 
tracks (see Section IIII A[) with momentum in the range 
1.1 GeV/c < p < 1.5 GeV/c and examining the timing 
difference between the measured fiight-time and the time 
which is expected under the assumption that the parti- 
cles are pions. The resulting time distribution is shown 
in Fig. [21 Since the pions dominate the total yield in 
this momentum region, a narrow peak centered around 
t — ioxpcctcd ~ is observed. The other two broad peaks 
in Fig. [2] correspond to kaons and protons. A Gaussian 
distribution is fit to the pion peak and yields a resolution 
of ctbbc-tofw = 84 ± 1 ps for the BBC-TOFw system. 



x10' 
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FIG. 2; (color online) Timing difference T — T^, the differ- 
ence between the measured time in the TOFw and the time 
calculated assuming each candidate track is a pion. 

The excellent timing resolution allows for 4cr separa- 
tion in mass-squared reaching up to = 2.5 GeV/c for 
tt/K, and up to pt = 4 GeV/c for K/p. The PID is 
further extended in px by use of asymmetric cuts around 
the centroids of the mass-squared distributions. 



TABLE I: The main characteristic parameters of TOFw and 
AGO 



TOFw 



ACC 



Arj 
A(j> (rad) 

Radial distance (cm) 
number of cells 
cell size (cm^) 



(-0.35, 0.35) 
(-0.061, 0.110) 
(0.503, 0.674) 
481.36 

512 
37 X 2.8 



(-0.35, 0.35) 
(-0.108, 0.156) 

449.4 
160 
11.95 X 23.10 



0.3 cm. The readout is double-sided, which allows The ACC is used in conjunction with the TOFw to 
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aid the PID at high px- It is situated in the West spec- 
trometer arm in front of the TOFw detector. The ACC 
is a Cerenkov radiation detector with a relatively high 
index of refraction (n = 1.0113, jth = 8.5), which means 
that light is produced at relatively low momenta. The 
threshold for radiation is 1.0 GeV/c for pions, 3.0 GeV/c 
for kaons, and 6.0 GeV/c for protons. The combined 
ACC-TOFw information allows for tt/K separation up to 
Pt = 4 GeV/c, and K/p separation up to px = Q GeV/c. 
The main characteristic parameters of TOFw and ACC 
can be found in the Table HI 

III. ANALYSIS METHOD 
A. Event and Track Selection 

The results reported here are obtained from an anal- 
ysis of 4.8 X 10^ minimum bias events obtained during 
the 2007 running period. The minimum-bias trigger is 
defined by a coincidence between North and South BBC 
signals and an energy threshold of one neutron in both 
the North and South zero-degree calorimeters The 
collision vertex z is constrained to \z\ < 30 cm of the 
origin of the coordinate system. 

Charged tracks are selected based on the track quality 
information from the tracking system (DC-PCl). The 
tracks are then projected to the outer detectors and con- 
firmed by requiring that the closest hit to the track 
projection is within certain spatial windows in cj) and 
z. The distributions for the distance between the clos- 
est hit and projection in the azimuthal and z directions 
are fitted with a double Gaussian function, one Gaus- 
sian function is for the signal distribution and the other 
for the background. For < 3 GeV/c, hits are re- 
quired to match the TOFw and the PC3 to within 2a 
from the signal's Gaussian distribution in 4> and z. For 
Pt > 3 GeV/c, hits are required to match the PC2 and 
the PC3 to within 3cr and the TOFw to within 2a in 
and z. Background from 7-conversions is further reduced 
by applying a RICH veto. For the pions, this veto only 
works for px < 5 GeV/c since pions with px higher than 
that will fire the RICH. To evaluate the residual back- 
ground, remaining after these selections, the background- 
to-signal (B/S) ratios from the double Gaussian function 
fitting within the samples selected for the analysis are 
examined. For px < 3 GeV/c the background com- 
prises less than 1% of the selected tracks. At higher 
px the background increases, reaching B/S sa 7% for 
5.5 GeV/c < Px < 6.0 GeV/c in the 0-20% centrality 
bin. 



B. Particle Identification 

The particles are identified by their mass, based on 
measurements of the momentum, the time-of-fiight to 
the TOFw detector, and the path-length along the 



trajectory. PID selections are performed by applying 
momentum-dependent cuts in mass-squared. The mass- 
squared distributions are fit with a 3-Gaussian function 
corresponding to pions, kaons, and protons. The cor- 
responding widths and centroids are extracted from the 
data as a function of transverse momentum. In the cal- 
ibration process, we ensure that the centroids of these 
distributions do not move as a function of px and that 
the widths vary as expected from the known momen- 
tum and timing resolution of the detector. We then se- 
lect a sample from each particle species aiming for at 
least 90% purity in PID. The high purity of the sam- 
ple will allow us to measure the V2 of selected par- 
ticles accurately and minimize the uncertainty result- 
ing from PID contamination. At lower transverse mo- 
menta {px < 2.5 GeV/c), the 2a bands centered around 
each particle's do not overlap, thus symmetric cuts, 
ml — 2a < m? < 111^ + 2a, allow for PID with high pu- 
rity. In the range 2.5 GeV/c < pt < 3 GeV/c, the t:/K 
separation is achieved by excluding the particles that lie 
within 2(7 of the centroid of the mass-squared distribu- 
tion of another particle. This procedure is demonstrated 
in Fig. m where the PID selections for tt, K, and p are 
shown with the hatched areas in the plot. The Gaussian 
fits to the individual peaks (dashed-line curves) and 
the combined fit to the entire distribution (solid line) 
are also shown. 




m2 (GeV^/c") 

FIG. 3: (color online) The mass-squared distribution 
measured by TOFw in the pr region 2.9 GeV/c < 
Pt < 3.0 GeV/c. The hatched areas show the pion, kaon, and 
proton selections, from left to right. The dashed lines show 
Gaussian fits to the individual peaks, while the solid line 
represents a combined fit to the distribution including the 
pions and kaons. 

At higher transverse momentum 3 GeV/c < 
Pt < 6 GeV/c, the lower range of the pion distribu- 
tion remains relatively unaffected by contamination from 
kaons and protons. Therefore, a pion sample with purity 
better than 90% can be selected based on information 
from the TOFw alone, by applying the cuts indicated 
in panels (a) and (c) of Fig. 21 and listed in Table HIl 
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FIG. 4: (color online) The mass-squared distribution in the TOFw without (left panels) and with (right panels) the ACC photon 
yield (accph) cuts for different pr regions. The hatched areas show the cuts used for pion, kaon, and proton selections. The 
distribution is fit with a 3-Gaussian function (solid line). The individual Gaussian distributions corresponding to vr, K, and p 
and are as dashed lines. 
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FIG. 5: (color online) The mass-squared distribution in the TOFw without (left panels) and with (right panels) the ACC 
photon yield (accph) cuts for different pr regions. The distribution is fit with the sum of three empirical distribution 
functions that are propagated from sampling a Landau shape momentum distribution as described in the text (solid line). The 
individual distributions corresponding to tt, K, and p are shown with dashed lines. 
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For kaon and proton identification at pt > 3 GeV/c, 
the ACC is used in conjunction with the TOFw detec- 
tor, as shown in Fig. |4j panels (b) and (d). The turn- 
on momenta of the ACC for pions, kaons, and protons 
are 1.0 GeV/c, 3.0 GeV/c, and 6.0 GeV/c, respectively. 
This turn-on is gradual, with the number of photons reg- 
istered per photomultiplier tube (PMT) growing up to 15 
for pions, and 10 (kaons and protons) as the hadrons ex- 
ceed their respective threshold momentum by « 1 GeV/c. 
With this information, the photon yield from the ACC 
can be used as a rejection veto based on whether it is 
"on" (accph > 5) or "off" (accph < 5). Due to the occu- 
pancy effects in the ACC as well as the spatial resolution 
of track projection to the ACC, the pions cannot be re- 
jected completely. The effect of this veto cut is demon- 
strated in Fig. 01 panels (b) and (d). The pion rejection 
by the ACC in combination with asymmetric rri^ cuts, 
which are indicated here and listed in Table HIl allow for 
kaon and proton PID up to px of 4 and 6 GeV/c, respec- 
tively. 

We use the Gaussian fits to the mass-squared dis- 
tributions to estimate the PID purity in the selected 
771^ regions. This is straightforward at pt < 3 GeV/c, 
where the peaks associated with each particle are well 
defined. At higher px the uncertainties are larger, since 
the pion and kaon peaks merge and the individual yields 
are not well constrained. We have checked that for 
Pt > 4 GeV/c, after efficiency corrections, the K/ir ra- 
tio obtained from our fits is consistent with the mea- 
surements of the K/tt {Kg/ir) ratio by the STAR ex- 
periment within the statistic and systematic uncertain- 
ties. At Pt = 5.22 GeV/c, the K/tt ratio is reported as 
0.326 ± 0.013(stat) ± 0.134(syst) and the Xy-K ratio 
is reported as 0.435 ± 0.022(stat) ± 0.072(syst) in p+p 
collisions at y^SNN — 200 GeV; the ratios in p+p and 
Au-|-Au collisions are similar |44|-|46| . In our study, the 
kaon contamination in the pion sample is relatively in- 
sensitive to the kaon yield. For example, if we artificially 
increase the kaon yield by 30%, the contamination in the 
pion sample increases from 7% to 9%. 

The distributions are not strictly Gaussian shape, 
but have tails extending to the higher mass region. This 
effect is not noticeable at low pt but comes into promi- 
nence at intermediate and high pT- Hadrons coming 
from resonance decays may survive the tracking cuts but 
will have misreconstructed momentum and contribute to 
this high mass tail. The total momentum distribution of 
hadrons, including those from resonance decays, is much 
closer to a Landau distribution than a Gaussian distri- 
bution. To get an estimate of the possible PID contami- 
nation in this case, we have fit the distribution with 
an empirical function that was determined by sampling 
a momentum distribution with a Landau shape instead 
of a Gaussian. This empirical distribution is found to 
give a much better fit than a simple 3-Gaussian function 
and it gives a good description of the high mass tails. 
Finally, we reevaluate the PID contamination with this 
empirical function. An example of these fits is shown in 



Fig. [5j The tail of distribution is well described by 
the empirical pion, kaon, and proton functions which 
are presented with different dashed lines. In this case, at 
high Pt the contamination from kaons and pions in the 
proton sample increases to 9% from 1% in the case of the 
Gaussian fits. 

The PID purity for each particle species estimated in 
different pT ranges is listed in Table HT] These estimates 
reflect the values obtained for the 0-20% central Au-l-Au 
collisions and are meant to provide lower limits for the 
measurements presented here. The purity in more pe- 
ripheral collisions was found to be slightly better. 



C. Measurement of V2 

The measurement of the anisotropy parameter V2 aims 
to determine the event- by- event particle azimuthal cor- 
relation with the reaction plane of the collision. The 
true reaction plane, which is defined as the plane formed 
by the impact parameter b and the beam direction, is 
not known experimentally. In addition, there exist other 
sources of correlations in azimuth, such as the correla- 
tions from resonance decays, jets, and quantum effects. 
These correlations, which are not related to the reac- 
tion plane, are called nonflow correlations. The goal is 
to determine the second coefficient in the Fourier expan- 
sion V2 of the particle azimuthal distribution with respect 
to the reaction plane with minimal effects from nonflow 
correlations. To estimate the reaction plane angle ^'/jp, 
we employ the event-plane method (40| , in which the sec- 
ond harmonic azimuthal anisotropy signal determines the 
event-plane angle ^ based on hits registered in one of the 
event-plane detectors: RxNP or MPC. For an ideal de- 
tector, the measured distribution of event-plane angles 
should be isotropic. However, the actual measurement 
is usually affected by flnite acceptance and nonuniform 
efflciencies. We apply a standard event-plane flattening 
technique [H, [T9l,l40l |47| to remove the residual nonuni- 
formities in the distribution of event-plane angles. The 
accuracy with which the event-plane angle can be de- 
termined depends on the strength of the V2 signal and 
the multiplicity of the events in each centrality class. It 
is maximal for midcentral events, where both of these 
quantities are relatively large. The V2{pt) measurement 
is performed by correlating the particle azimuthal angle tp 
with the second harmonic event-plane angle 5*, and cor- 
recting the observed signal for the event-plane resolution 
as follows: 



V2 



(cos(2(^ - M>))) 
Res(*) 



(2) 



Here the brackets () indicate an average over all particles 
in all events. 

Since the true reaction plane angle is not directly mea- 
surable, the resolution correction is estimated using sub- 
event techniques [i^l • There are several different options 
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TABLE II: The particle identification cuts in TOFw and ACC witii PID purity in Au+Au collisions for the centralities 0-20% 
and 20-60%. 



Particle 


Pt range 


iUrw Cut 


AOO Out 


Purity 












U— ZU /o 


zU— DU /o 


pion 








IN one 


QQ% 


QQ% 






veto on 771 dz 1(7^2 










[3.0, 5.0) 


[—1.0, 0.0] 


None 


95% 


96% 




[O.U, D.Oj 


r 1 A All 
[-1.0, -O.IJ 


None 


yi/o 


9270 


kaon 


< 6 




None 


9o7o 


99% 






veto on ± 2cr„2 










[3.0, 3.5) 


[0.1,0.5] 


accph < 5.0 


94% 


95% 




[3.5,4.0) 


[0.1,0.5] 


accph < 5.0 


91% 


92% 


proton 


< 3 


ml ± 2cr„2 


None 


99% 


99% 




[3.0,4.0) 


[0.6, 1.3] 


accph < 5.0 


97% 


98% 




[4.0,5.0) 


[0.7, 1.3] 


accph < 5.0 


95% 


96% 




[5.0, 6.0) 


[0.7, 1.7] 


accph < 5.0 


91% 


92% 



in using the sub-event techniques. The present analysis 
uses the two sub-event and the three sub-event methods. 
These methods are compared to evaluate the systematic 
uncertainties associated with the event-plane resolution 
corrections. 
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FIG. 6: (color online) Panel (a) shows the event-plane reso- 
lution as a function of centrality for the RxNP and the MPC 
detectors. Panel (b) shows the ratio of the event-plane reso- 
lution obtained from two sub-events and three sub-events as 
a function of centrality. 

The RxNP and the MPC detectors each have two sub- 
detectors, North and South, which are positioned sym- 
metrically around the origin of the nominal collision point 
with equal acceptance in pseudorapidity. Thus, they pro- 



vide a natural two sub-event division. The correlation be- 
tween the event-plane angles determined from the North 
and South sub-detectors, n and s^ allows for the es- 
timate of the resolutions corrections as follows: 



Res(*Ar) = Res(^'s) = V(cos2 (^-g - *iv)) (3) 

The resolution correction can also be expressed analyti- 
cally [13 as: 



Res(*) = (cos2 - -^Rp)) 



(4) 



where /q and 1\ are modified Bessel functions. The pa- 
rameter X — W2 V2M, where M is the number of particles 
used to determine the event-plane, describes the disper- 
sion of the flow vector. With the use of Equation [3] and 
Equation m we obtain the sub-event parameters xs and 
Xn- Subsequently, to optimize the event-plane resolu- 
tion, the two sub-events are combined, and the full event 
parameter is taken as x = V^XS = ^J^Xn- This pro- 
cedure relies on the two sub-events being equal in mul- 
tiplicity and registering the same size W2 signal, which 
may not be the case experimentally. To avoid this un- 
certainty, we also use a three sub-events technique to de- 
termine the event-plane resolution with Equation [S] 140|] . 
To determine the event-plane resolution of RxNP detec- 
tor (sub-event A), we employ information from the North 
and South portions of the MPC detector (sub-events B 
and C). In turn, to estimate the resolution of the MPC 
detector, the North and South portions of the RxNP de- 
tector are used to provide sub-events B and C. 



Res(*^ 



(cos2(«'A - *i?p)) 



10 



0.30 
0.25 
0.20 

CM 

> 0.15 
0.10 
0.05 



O 1-15 
Q_ 1.10 
^ 1.05 
1.00 

^ 0-95 
P 0.90 
^ 0.85 



. . . 1 . . . . 1 . 1 . . 1 . . . . 1 . 1 . . 1 . . . . 1 . I_ 

r (a) pion 0-60% 

7 ORXN -. 
'-_ •MPC J 

- T - 

:i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 h 1 1 1 1 1 1 1 1 1: 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

r (b) kaon 
• 

j 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 M ! 


r (c) proton " 

: 

•• 
• 
• 

j 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1. 1 1 1 1 1 1 ! 


I „ ! 


f(e) 


E"* , t ! 


[ AMUMf^^ ^ j ^ 


i iAii^V^, 

;. ... 1 .... 1 .... 1 .... 1 .... 1 .... 1 .. ; 


1 1 



2 3 4 5 6 
(GeV/c) 



1 2 3 4 5 6 
(GeV/c) 



2 3 4 5 6 
p^ (GeV/c) 



FIG. 7: (color online) The upper panels show the azimuthal anisotropy V2{pt) of pions (a), kaons (b), and protons (c) in 
the 0-60% centrality class measured with respect to event-planes determined by the MPC (closed symbols) or the RxNP 
(open symbols) detectors. The event-plane resolution is estimated by the three sub-events method. The ratio of (RxNP) to 
U2(MPC) is shown in the lower panels as a function of pT for pions (d), kaons (e), and protons (f). 



/ (cos 2 {^A - ^b)) (cos 2 - ^c)) , . 

y (cos2(*B-*c)> 

The event-plane resolution for the RxNP (circles) and 
the MPC (triangles) detectors obtained with the above 
procedures are shown as a function of the event centrality 
in Fig. El^a). The results show the expected trend, with 
maximal resolution for the 20-30% centrality class where 
both the event multiplicity and the V2 signal are large, 
and a decrease for the more central events (due to lower 
V2 strength) , and more peripheral events (due to smaller 
multiplicity) . Figure IHJb) shows the ratio of the results 
obtained with the two sub-event and the three sub-event 
techniques. The results for the RxNP detector (closed 
symbols) agree to within 2%. A larger difference (up to 
4%) is observed for the MPC detector (open symbols) 
which is mainly due to the asymmetric pseudorapidity 
coverage of the MPC. The event-plane resolution from 
three sub-events method is used to correct the V2 mea- 
surement. 

From Fig. [5] it is evident that the RxNP detector has 
better resolution for the event-plane angle, as well as 
smaller systematic uncertainty in the event-plane deter- 
mination than the MPC detector. Therefore, it is de- 
sirable to use the RxNP for the V2 measurement. One 
possible disadvantage of the RxNP over the MPC de- 
tector is the smaller pseudorapidity separation from the 
central spectrometer (|?7| < 0.35), which makes the V2 
measurement more susceptible to nonflow correlations 
caused by jets. Since the results presented here aim to 



study the high pT azimuthal anisotropies of identified 
charged hadrons and the limits of Ug scaling, it is partic- 
ularly important to minimize such effects. To evaluate 
the nonflow contributions we examine the V2(j>t) distri- 
butions for pions, kaons, and protons measured using the 
MPC and the RxNP detectors independently. For the 
Pt < Q GeV/c, a previous study indicated that nonflow 
effects are small for the event-plane measured by the BBC 
detectors, which have a pseudorapidity coverage similar 
to that of the MPC detectors [4^. Figure [7] shows the 
results in the 0-60% centrality range for each particle 
species (upper panels), and the ratio of the results ob- 
tained with the two event-plane detectors (lower panels) . 
Non-flow correlations are expected to enhance the mea- 
sured V2 signal for the detector which is more affected, 
especially in the higher pT range. We do not find any ev- 
idence for a significant increase in nonflow contributions 
in the measurement based on the RxNP detector. 

Based on these considerations, the results presented 
in Section IIVI are based on the reaction plane measured 
solely by the RxNP, taking advantage of its better event- 
plane resolution in comparison to the MPC. 



D. Systematic Uncertainties in V2 

The systematic uncertainties in the V2 measurement 
obtained with the RxNP detector can be broadly charac- 
terized according to the following categories: 1) event- 
plane resolution corrections; 2) event-plane measured 
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TABLE III: Systematic uncertainties given in percent on the V2 measurements. 



Error Sources 



Species 



- 20% 



20 - 60% 



Type 



Event-plane resolution 
Event-plane detectors 

Background 



PID 



Acceptance 
and run-by-run 



pion 

kaon 
proton 

pion 

kaon 
proton 



2% 

3% in pT 1-3 
5% in PT 3-6 
1% in pt 1-4 
4% in PT 4-6 
1% in PT 1-4 
1% in pt 1-4 
5% in PT 4-6 



3% in PT 3-4 
5% in PT 4-6 



GeV/c 
GeV/c 
GeV/c 
GeV/c 
GeV/c 
GeV/c 
GeV/c 

negligible in 
2% in Pt 

negligible in 
2% in Pt 

negligible in 
GeV/c 
GeV/c 



3% in 
5% in 
1% in 
3% in 
1% in 
1% in 
3% in 
PT 1-3 GeV/c 
3-6 GeV/c 
PT l-3GeV/c 
3-4 GeV/c 
PT 1-3 GeV/c 
2% in 
3% in 



2% 
PT 1-5 
PT 5-6 
Pt 1-4 
Pt 4-6 
Pt 1-4 
Pt 1-4 
Pt 4-6 



GeV/c 
GeV/c 
GeV/c 
GeV/c 
GeV/c 
GeV/c 
GeV/c 



Pt 3-4 
Pt 4-6 
3% 



GeV/c 
GeV/c 



C 
B 

A 

A 
A 

A 

A 

A 



from different detectors; 3) V2 from background tracks; 
4) PID purity; and 5) acceptance and run-by-run depen- 
dencies. 

The uncertainties stemming from the event-plane reso- 
lution corrections are independent of particle species and 
Pt- They are found to be around 2% for all centralities 
by studying the event-plane resolution difference for the 
RxNP with the two and three sub-event methods. 

The uncertainties from event-planes measured with dif- 
ferent detectors (RxNP, MPC) are found to be indepen- 
dent of the particle species, by comparing the results from 
RxNP and MPC. In the 0-20% centrality class, we assign 
a 3% systematic uncertainty for px < 3 GeV/c and a 5% 
systematic uncertainty for px > 3 GeV/c. In the 20-60% 
centrality class, we assign a 3% systematic uncertainty 
for Pt < 5 GeV/c and a 5% systematic uncertainty for 
Pt > 5 GeV/c. 

Background tracks that are not removed by the track- 
ing and PID selections outlined in Sections llII Al and llll Bl 
may influence the measured V2 if they carry a signal 
which is different from the particle of interest. The back- 
ground tracks may come from decays, 7-conversions, or 
false track reconstruction. The backgrounds are central- 
ity dependent, and may also have Pt and hadron species 
dependence. A sample of background-dominated tracks 
was selected based on the normalized distance between 
the hits registered in the TOFw detector and the track 
projections. Specifically, a Aa-lOa window in the z- 
direction was utilized. The azimuthal anisotropy of the 
background was then measured following the procedure 
used for the signal. For pT ^ 3 GeV/c, the V2 of the back- 
ground is similar to that of the pion, but it decreases at 
higher pt down to about 60% of the V2 of the pion (or 
30% of the V2 of the proton) for p-r ~ 6 GeV/c in the 0- 



20% centrality class. For pions, the resulting systematic 
uncertainties in V2 are of the order 1% for < 4 GeV/c 
and reach up to 4% (3%) for pT ~ 6 GeV/c for centrality 
0-20% (20-60%). For protons, the resulting systematic 
uncertainties in V2 are of the order 1% for pT < i GeV/c 
and reach up to 5% (3%) for pT ~ 6 GeV/c for centrality 
0-20% (20-60%). 

The systematic uncertainties in V2 resulting from 
hadron misidentification are based on the PID purity esti- 
mates listed in Table [Til and the size of V2 of each species. 
For example, at pT — 6.0 GeV/c for 0-20% centrality, 
the protons purity is around 91% and the V2 of pions and 
kaons are around 50% of that of the proton. We assign 
a 5% systematic uncertainty attributable to this effect. 
For Pt < 3 GeV/c, the uncertainties in V2 due to PID 
contamination are negligible for all particle species. At 
higher pT the uncertainties in V2 remain below 2% for 
kaons and pions; for protons with pT > 4: GeV/c these 
uncertainties reach up to w 5% (3%) for centrality 0-20% 
(20-60%). 

Additional systematic checks were performed using dif- 
ferent subsets of the detector, and data obtained with dif- 
ferent magnetic field configurations. Differences of order 
8% (3%) were found for the 0-20% (20-60%) centrality, 
which have weak pT and particle species dependence. 

Table [nil lists the summary of all these systematic un- 
certainties which are categorized by the types: 

(A) point-to-point error uncorrelated between pT bins, 

(B) Pt correlated, all points move in the same direction 
but not by the same factor, 

(C) an overall normalization error in which all points 
move by the same factor independent of pt- 
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FIG. 8: (color online) Identified hadron V2 in central (0-20% centrality, left panels) and midcentral (20-60%, right panels) 
Au+Au collisions at ^snn = 200 GeV. Panels (a) and (b) show V2 as a function of transverse momentum pr- Panels (c) and 
(d) show the quark- number-scaled 112 {v2/nq) as a function of the kinetic energy per quark, KEr/nq. Panels (e) and (f) show 
v^/uq as a function of transverse momentum per quark, pT/riq. The V2 of all species for centrality 0-20% has been scaled up by 
a factor of 1.6 for better comparison with results of 20-60% centrality. The error bars (shaded boxes) represent the statistical 
(systematic) uncertainties. The systematic uncertainties shown are type A and B only. 



IV. RESULTS AND DISCUSSION 



The results for V2 of identified pions, kaons, and pro- 
tons are presented in Fig. [51 the results in central col- 
lisions (0-20%) are presented in panels (a), (c), and (e) 
and the results in noncentral collisions (20-60%) are pre- 
sented in panels (b), (d), and (f). The symbols repre- 
senting the different particle species are closed triangles 
for pions, open squares for kaons, and closed circles for 



protons. In order to better compare between two cen- 
tralities, the V2 of all species in the 0-20% centrality has 
been scaled up by a factor of 1.6. The error bars (shaded 
boxes) represent the statistical (systematic) uncertain- 
ties. The systematic uncertainties shown are type A and 
B only. Not shown are the type C systematic uncertain- 
ties, which are from the event-plane resolution, geomet- 
rical acceptance, and run-by-run dependence are around 
8.5% (3.5%) for 0-20% (20-60%) centrality for all species 
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at all values of pt- 

Panels (a) and (b) of Fig. [5] show V2{pt)- For both 
centrality selections, the V2 values of pions and kaons are 
very similar in intermediate pt range (2-4 GeV/c) , where 
the measured V2 is maximal and is relatively independent 
of transverse momentum. Above pt ~ 4 GeV/c the pion 
V2 gradually decreases to a value which is comparable to 
the signal measured at ~ 1 GeV/c. In contrast, the 
proton V2{pt) has a shape which is centrality dependent. 
In central collisions (0-20%) the proton V2 rises up to 
Pt ~ 3.5 GeV/c and then saturates at a value higher than 
the V2 of pions. For noncentral collisions, the behavior 
is different: a decrease is observed in the proton V2 for 
Pt > 4 GeV/c leading to near equal V2 signals for pions 
and protons at ~ 6 GeV/c. 

The use of the KE-r variable was introduced in Refer- 
ence which is found to better represent the number 
of quark scaling behavior than pT at lower px- In panels 
(c) and (d) of Fig. [8] the V2 signals have been scaled by 
the number of constituent quarks Ug in the hadrons and 
are shown as a function of the transverse kinetic energy 
per quark KEx/riq. A very different behavior is observed 
in central (Fig. [Hljc)) and in noncentral (Fig. [Sjd)) col- 
lisions. In the measured pt range, a universal behavior 
is seen in the central collisions within the statistical and 
systematic uncertainties, but not in the noncentral colli- 
sions, where the V2/ng of protons falls below that of the 
mesons for KEr/ng > 1 GeV. This is the range where 
the proton V2{pt) begins falling in noncentral collisions 
but remains relatively constant in central collisions. 

On the other hand, it is widely accepted that the 
relevant scaling variable for quark-recombination is the 
transverse momentum per quark, since it is the momen- 
tum and not the energy that is additive in the recom- 
bination models. Therefore, to examine the Uq scaling 
in the recombination regime we show the quark-number- 
scaled V2 as a function of pr/nq in panels (e) and (f) 
of Fig. m For central collisions (Fig. EJe)), the univer- 
sal behavior appears to remains valid within the statis- 
tical and systematic uncertainties. Since the changes in 
V2 are relatively small at higher pr, shifting the cc-axis 
from KEr/uq to pr/nq does not change the shape of the 
curves significantly. For noncentral collisions (Fig.lSff)), 
the proton data are systematically below the pion data 
at all pr/uq, although they are at the edge of the system- 
atic uncertainties for pT/nq < 1.3 GeV/c, which corre- 
sponds to KEx/riq < 1 GeV/c. We note that despite this 
systematic offset, the Uq scaling makes the shape of the 
pion and proton curves very similar below the breaking 
point. Above that point, quark recombination is clearly 
violated. 

Some model calculations \3d\ have shown that the 
breaking of Uq scaling occurs at the transition between 
purely thermal and thermal-|-shower recombination. In 
the 50-60% centrality class this can happen for values 
of KEt as low as KFix/nq « 0.5 GeV, while in the 
0-5% centrality class this occurs at values as high as 
KEr/n-g « 1.5 GeV. Similar features have been observed 



in the data presented in this paper. On the other hand, 
for pions and protons, the nuclear modification factors 
(Raa), which are used to quantify the amount of par- 
tonic energy loss in the medium, have been found to be 
consistent with each other for pT > 5 GeV/c [46l l49l - [5]| . 
This indicates that a simple interplay between recombi- 
nation and jet energy loss is not enough to explain the 
V2 and Raa of pions and protons in Au-|-Au collisions 
in this Pt region. Additional considerations may include 
the nonAbelian nature of jet energy loss js^], the quark 
versu s gl uon fragmentation production of pions and pro- 
tons |53l - l55| . and jet chemistry effects such as enhanced 
parton splitting [sl] and jet conversion [s^l- Detailed 
model calculations that take all of these effects into ac- 
count are not yet available, and it is an open question 
whether doing so is enough for an adequate interpreta- 
tion of the Pt V2 and -Raa of pions and protons. 



To further investigate the centrality dependence of the 
Uq scaling breaking, results with finer centrality bins are 
shown in Fig. |9l The quark number scaled V2 {v2/nq) of 
pions, kaons, and protons are shown as a function of the 
kinetic energy per quark KE^/n^ in 0-10% (panel (a)), 
10-20% (panel (b)), 20-40% (panel (c)), and 40-60% cen- 
trality (panel (d)). The error bars (shaded boxes) repre- 
sent the statistical (systematic) uncertainties. The sys- 
tematic uncertainties shown are type A and B only. Not 
shown are the type C systematic uncertainties, which are 
from the event-plane resolution, geometrical acceptance, 
and run-by-run dependence are around 10.5% (3.5%) for 
0-10% (40-60%). These resuhs with finer centrality bins 
show that the breaking of nq scaling has a clear centrality 
dependence. 



We also compare our results with the existing V2 results 
for Kg and A as measured by the STAR collaboration us- 
ing the event-plane method [2l| in the 0-10% and 10-40% 
centrality classes, which are shown in panel (a) and panel 
(b) of Fig. [TUl respectively. Since the event-plane and 
particles are measured in the same rapidity gap by the 
STAR detector in their event-plane method, the V2 values 
from STAR measurements are expected to be systemat- 
ically larger than those measured by PHENIX [H \^ 
due to nonflow effects. In the 0-10% centrality class, the 
V2 of pions and protons in this study are systematically 
lower than the V2 of Kg and A by 17% independent oipT, 
but they are within the systematic uncertainties. The Uq 
scaling appears to hold in this centrality class for each 
particle species. In the 10-40% centrality class, the V2 
of pions and protons are consistent with that of Kg and 
A in the overlapping KE^ region. While the presence of 
the scaling breaking is not clear in the Kg and A results, 
the improved precision and extended KEy reach of the 
present study unambiguously demonstrates the breaking 
of Uq scaling in this centrality class. 
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FIG. 9: (color online) The quark-number-scaled V2 {v2/nq) of identified hadrons are sliown as a function of the kinetic energy 
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FIG. 10: (color online) The quark-number-scaled V2 (v^lfiq) of identified hadrons are shown as a function of the kinetic energy 
per quark, KEr/ng in 0-10% centrality (panel (a)) and 10-40% centrality (panel (b)) in Au-f Au collisions at -yijvjv = 200 GeV. 
The V2 of A and K% are measured by STAR collaboration [SJ] . The error bars (open boxes) represent the statistical (systematic) 
uncertainties. The systematic uncertainties shown on the results from this study are type A and B only. 
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V. SUMMARY 

We have presented a high-statistics study of baryon 
and meson azimuthal anisotropy V2 measured up to px of 
6 GeV/c as a function of centrahty in ^/snn = 200 GeV 
Au+Au collisions. The scaling is found to exhibit 
strong dependence on the collision centrality. Signifi- 
cant deviations from Uq scaling are found in noncentral 
collisions, starting from the 10-20% centrality class, as 
KEt/hq > 0.7 GeV. These results indicate that parton 
fragmentation and the associated energy loss may play an 
important role in generating the azimuthal anisotropy of 
particle emission. Conversely, in central collisions, such 
as the 0-10% centrality class, the universal Uq scaling 
appears to hold to KE^/ng — 1.5 GeV, supporting par- 
ton recombination as the dominant mode of particle pro- 
duction at intermediate transverse momentum in central 
Au+Au collisions at top RHIC energy. 

ACKNOWLEDGMENTS 

We thank the staff of the Collider-Accelerator and 
Physics Departments at Brookhaven National Labora- 
tory and the staff of the other PHENIX participating 
institutions for their vital contribiitions. We acknowl- 
edge support from the Office of Nuclear Physics in the 
Office of Science of the Department of Energy, the Na- 
tional Science Foundation, Abilene Christian University 



Research Council, Research Foundation of SUNY, and 
Dean of the College of Arts and Sciences, Vandcrbilt Uni- 
versity (U.S. A), Ministry of Education, Cultiue, Sports, 
Science, and Technology and the Japan Society for the 
Promotion of Science (Japan), Conselho Nacional de De- 
senvolvimento Cientifico e Tecnologico and Fundagao de 
Amparo a Pesquisa do Estado de Sao Paulo (Brazil), Nat- 
ural Science Foundation of China (P. R. China), Min- 
istry of Education, Youth and Sports (Czech Repub- 
lic), Centre National de la Recherche Scientifique, Com- 
missariat a I'Energie Atomique, and Institut National 
dc Physique Nucleaire ct dc Physique des Particulcs 
(France), Ministry of Industry, Science and Tekhnolo- 
gies, Bundesministerium fiir Bildung und Forschung, 
Deutschcr Akademischcr Austausch Dienst, and Alexan- 
der von Humboldt Stiftung (Germany), Hungarian Na- 
tional Science Fund, OTKA (Hungary), Department of 
Atomic Energy and Department of Science and Technol- 
ogy (India), Israel Science Foundation (Israel), National 
Research Foundation and WCU program of the Ministry 
Education Science and Technology (Korea), Ministry of 
Education and Science, Russian Academy of Sciences, 
Federal Agency of Atomic Energy (Russia), VR and the 
Wallenberg Foundation (Sweden), the U.S. Civilian Re- 
search and Development Foundation for the Independent 
States of the Former Soviet Union, the US-Hungarian 
Fulbright Foundation for Educational Exchange, and the 
US-Israel Binational Science Foundation. 



[1] K. Adcox ct al. (PHENIX Collaboration), Nucl. Phys. A 

757, 184 (2005). 
[2] J. Adams et al. (STAR Collaboration), Nucl. Phys. A 

757, 102 (2005). 
[3] B. B. Back et al. (PHOBOS Collaboration), Nucl. Phys. 

A 757, 28 (2005). 
[4] I. Arscnc ct al. (BRAHMS Collaboration), Nucl. Phys. 

A 757, 1 (2005). 
[5] E. V. Shuryak, Nucl. Phys. A 750, 64 (2005). 
[6] M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 

(2005). 

[7] U. W. Heinz, arXiv:0901.4355 (2009). 
[8] P. Romatschkc, Int. J. Mod. Phys. E 19, 1 (2010). 
[9] D. A. Teaney, arXiv:0905.2433 (2009). 
[10] P. K. Kovtun, D. T. Son, and A. O. Starinets, Phys. Rev. 

Lett. 94, 111601 (2005). 
[11] S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 

C 76, 034904 (2007). 
[12] S. Afanasicv ct al. (PHENIX Collaboration), Phys. Rev. 

C 80, 054907 (2009). 
[13] A. Adare et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 105, 142301 (2010). 
[14] S. A. Bass, C. Gale, A. Majumder, C. Nonaka, G.-Y. 

Qin, ct al., Phys. Rev. C 79, 024901 (2009). 
[15] X. N. Wang, Phys. Rev. C 63, 054902 (2001). 
[16] S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 91, 182301 (2003). 



[17] J. Adams ct al. (STAR Collaboration), Phys. Rev. Lett. 

92, 052302 (2004). 
[18] A. Adare et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 98, 162301 (2007). 
[19] S. Afanasiev et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 99, 052301 (2007). 
[20] B. Abolcv ct al. (the STAR Collaboration), Phys. Rev. 

C 75, 054906 (2007). 
[21] B. Abelev et al. (STAR Collaboration), Phys. Rev. C 77, 

054901 (2008). 

[22] R. C. Hwa and C. B. Yang, Phys. Rev. C 67, 034902 
(2003). 

[23] R. J. Fries, B. MuUer, C. Nonaka, and S. A. Bass, Phys. 

Rev. Lett. 90, 202303 (2003). 
[24] R. J. Fries, B. MuUer, C. Nonaka, and S. A. Bass, Phys. 

Rev. C 68, 044902 (2003). 
[25] V. Greco, C. M. Ko, and P. Levai, Phys. Rev. Lett. 90, 

202302 (2003). 

[26] D. Molnar and S. A. Voloshin, Phys. Rev. Lett. 91, 
092301 (2003). 

[27] S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 91, 172301 (2003). 
[28] S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 

C 69, 034909 (2004). 
[29] B. MuUer, R. J. Pries, and S. A. Bass, Phys. Lett. B 618, 

77 (2005). 

[30] C. B. Chiu, R. C. Hwa, and C. B. Yang, Phys. Rev. C 



16 



78, 044903 (2008). 
1] G. Ferini, M. Colonna, M. Di Toro, and V. Greco, Phys. [45 

Lett. B 670, 325 (2009). 
2] K. Dusling, G. D. Moore, and D. Teaney, Phys. Rev. C [46 

81, 034907 (2010). 
3] R. A. Lacey, A. Taranenko, R. Wei, N. Ajitanand, [47 

J. Alexander, et al., Phys. Rev. C 82, 034910 (2010). 
4] O. Fochler, Z. Xu, and C. Greiner, Phys. Rev. Lett. 102, [48 
202301 (2009). 

5] R. A. Lacey, N. Ajitanand, J. Alexander, X. Gong, J. Jia, [49 

et al., Phys. Rev. C 80, 051901 (2009). 
6] T. Hirano and Y. Nara, Phys. Rev. C 69, 034908 (2004). [50 
7] J. C. Dunlop, M. A. Lisa, and P. Sorensen, Phys. Rev. C 

84, 044914 (2011). [51 
8] K. Adcox et al. (PHENIX Collaboration), Nucl. Instrum. 

Meth. A 499, 469 (2003). [52 
9] E. Richardson et al. (PHENIX Collaboration), Nucl. In- [53 

strum. Meth. A 636, 99 (2011). 
0] A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58, [54 
1671 (1998). 

1] M. Chiu (PHENIX CoUaboration), AIP Conf. Proc. 915, [55 
539 (2007). 

[42] B. Bonner et al., Nucl. Instr. Meth. A 508, 181 (2003). [56 
[43] W. J. Llope et al., Nucl. Instrum. Meth. A 596, 430 

(2008). [57 
[44] J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 



92, 052302 (2004). 

J. Adams et al., (STAR CoUaboration), nucl-ex/0601042 
(2006). 

G. Agakishiev et al., (STAR Collaboration), 
larXiv: 1110.0579 (2011). 

S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 
Lett. 96, 032302 (2006). 

S. Afanasiev et al. (PHENIX CoUaboration), Phys. Rev. 
C 80, 024909 (2009). 

B. I. Abelev et al. (STAR Collaboration), Phys. Rev. 
Lett. 97, 152301 (2006). 

A. Adare et al. (PHENIX Collaboration), Phys. Rev. C 
83, 064903 (2011). 

R. Belmont (PHENIX Collaboration), Nucl. Phys. A 
830, 697c (2009). 

X. N. Wang, Phys. Rev. C 58, 2321 (1998). 

S. Albino, B. A. Kniehl, and G. Kramer, Nucl. Phys. B 

725, 181 (2005). 

S. Albino, B. A. Kniehl, and G. Kramer, Nucl. Phys. B 
803, 42 (2008). 

D. de Florian, R. Sassot, and M. Stratmann, Phys. Rev. 
D 76, 074033 (2007). 

S. Sapeta and U. A. Wiedemann, Eur. Phys. J. C55, 293 
(2008). 

W. Liu and R. J. Fries, Phys. Rev. C 77, 054902 (2008). 



